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Abstract. Using the method proposed in a previous paper, we calculate pulsar braking indices in the models 
with torque contributions from both inner and outer accelerating regions, assuming that the interaction between 
them is negligible. We suggest that it is likely that the inverse Compton scattering induced polar vacuum gap and 
the outer gap coexist in the pulsar magnetosphere. We include the new near threshold vacuum gap models with 
curvature-radiation and inverse Compton scattering induced cascades, respectively; and find that these models 
can well reproduce the measured values of the braking indices. 
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1. Introduction 

Why do pulsars spin down? Although it is generally ac- 
cepted that pulsars brake due to magnetodipole radiation, 
this simple model is questioned by the fact that the ob- 
served braking index n = ilCt/Ct^ {il is the angular veloc- 
ity of rotation) is smaller than 3. For the Vela pulsar, the 
index is n = 1.4, which is the smallest one observed. 

Actually, some efforts have appeared to find unusual 
mechanisms to understand the observed braking index. It 
was suggested previously that the observation oi n ^ 3 
may result from a force-free precession of a distorted neu- 
tron star (Macy 1974) or the existence of a companion 
star (Deminanski & Proszynski 1979). Peng et al. (1982) 
and Huang et al. (1982) suggested that neutrino and pho- 
ton radiation from superfluid neutron vortexes may con- 
tribute to pulsar spindown. However, this idea may have 
difficulty in interpreting n < 3 since this mechanism dom- 
inates for older pulsars with longer periods. Blandford & 
Romani (1988) interpreted the observed index by a multi- 
pole field and/or field evolution. Melatos (1997) accounted 
for the spin-down of three pulsars (the Crab, PSR B0540- 
69, and PSR B1509-58) in the light of a non-standard 
vacuum dipole model, where a "vacuum radius" is intro- 
duced phenomenologically. Recently, accretion torque has 
been suggested to explain the discrepancy (Morley 1993, 
Menou, Perna & Hernquist 2001). Alpar et al. (2002) pro- 
posed that both the usual magnetodipole radiation and 
the propeller torque applied by the debris disk formed 



Send offprint requests to: R.X. Xu, 
e-mail: rxxuObac . pku . edu . cn 



soon after supernova explosion should cause a pulsar to 
spin down (see also, e.g., Marsden et al. 2001), although 
no direct evidence exists to show that a normal pulsar has 
significant torque of this kind. 

An alternative effort, within the framework of stan- 
dard^ neutron stars and their magnetospheric emission 
models, was proposed by Xu & Qiao (2001; hereafter 
Paper I). They addressed that the magnetodipole radi- 
ation and the unipolar generator are two different sources 
powering a pulsar, and these two should be combined in 
studying pulsar spindown. They find that the calculated 
braking indices, although model-dependent, are less than 
3, and that a proposed model should be ruled out if it 
cannot result in an index as small as 1.4. 

In addition, the pulsar emission mechanism is still a 
great challenge even over 30 years after discovery. It is 
a common point that two classes of accelerators (i.e., 
the inner gaps and the outer gaps) may work on pulsar 
magnetosphere to reproduce its photon (radio to X-ray 
or 7— ray bands) emission, and two more subclasses of 
models, i.e., the space-charge-limited fiow models (e.g., 
Arons & Scharlemenn 1979; Harding & Muslimov 1998) 
and the vacuum gap models (e.g., Ruderman & Sutherland 
1975) for the inner accelerators appear in the literature. 
However, it is quite possible that both inner and outer 
gaps may coexist in a magnetosphere (Usov 2000; Paper 
I) since those two gaps may work in different field line 
regions in order to close the global electric current (e.g., 
Holloway 1975). 

In this paper, using the method developed in Paper I, 
we calculate the braking indices in the models where both 
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inner and outer accelerators lose the rotation energy, ne- 
glecting the possible interaction between those two kinds 
of gaps. Since a new polar gap model, the near threshold 
vacuum gap (NTVG; Gil & Mitra 2001) model, became 
available after Paper I was published, we will begin with 
computing the braking indices in this new model. It is 
found that the NTVG model can reproduce an index as 
small as 1.4, and this model is therefore possible to explain 
the 5 braking indices observed. However, if the inner and 
the outer gaps coexist, we suggest that the ICS-induced 
vacuum gap and the outer gap may work together at the 
same time, since the calculated braking indices in this case 
are more reasonable than that in other cases. 

2. Braking indices for the NTVG model 

The model of Ruderman & Sutherland (1975) is "user 
friendly" , but has two imperfect points: 1, the binding en- 
ergy problem of ions on the neutron star surface; 2, only 
half of the neutron stars are applicable (e.g., Xu, Qiao & 
Zhang 1999). Although these points can be overcome if 
some radio pulsars are bare strange stars (Xu et al 1999), 
Gil & Mitra (2001) argued that the RS-type vacuum gaps 
could also exist for neutron stars with multipolar surface 
magnetic fields (i.e., the actual surface strength is much 
higher than the dipolar surface component, and the radius 
of field curvature is much smaller than the neutron star 
radius) , since the ion cohesive energy becomes larger if the 
field is stronger. 

However, in a superstrong surface field, Bs > 0.1i?q ~ 
4.4 X lO'^^ G, the 7 photons should be converted into at 
or near the kinematic threshold (Gil & Mitra, 2001); the 
vacuum gap formed under this condition is thus called 
as the Near Threshold Vacuum Gap (NTVG) model. 
Gil & Melikidze (2002) present the gap heights and the 
gap potential drops for different pair-production mech- 
anisms (CR-NTVG: curvature-radiation induced NTVG; 
ICS-NTVG: inverse-Compton- scattering induced NTVG) 
and different estimates of the cohesive energies of surface 
iron ions (i.e., Abrahams & Shapiro 1991, hereafter AS91; 
and Jones 1986, hereafter J86). In this section, we calcu- 
late the braking index for the NTVG model in the frame 
proposed in Paper I. 

For pulsars with given pulsar radius R and polar mag- 
netic field at surface B, the braking index n can be written 
as 

nf) n dr] 

n = 3 + — = 3 H — -, (1) 

with 

77 ~ sin^ a + 5.4 X IQ-^R-^B^^^ cos^ an^'^Acj), 

where a is the inclination angle between the magnetic 
axis and the rotating axis. Re = i?/(10^cm), B12 = 
i3/(10^^G), and A<j) is the gap potential drop which is 
model dependent. 



For the curvature-induced NTVG model, the gap po- 
tential drop is (Gil & Melikidze 2002) 

A<^g^vG ^ 4 ^ l0\^/^pt/h'/'p-'/''P'_{f cgse, (2) 

with the parameter b in the form 

lCR _ oaO.08 7,0.57 ^0.32 p-l. 15 p-0. 5 
'^min — PS Pg ^ ^-15 ' 

where /3 = 52 for AS91 case and f3 = 1990 for J86 case, 
C ~ 0.85 is the typical value of the correction factor, k is 
the heat flow coefficient, P — 2tt/Q, is the pulsar period 
in seconds, and P_i5 = P/10~^^. For simplicity, one may 
take A: = 1.0 and — 1, and therefore has 

ryX"" - sin^ a + \?>.7p''-^^R^^-^^B^^ cos' a^-^-'^ (3) 

For the inverse-Compton-scattering-induced NTVG 
model, the gap potential drop is (Gil & Melikidze 2002) 

A(/.fcYG = 1.7xl0i°C°-''A:-0-iV6-''^"'^"'-''^'ri°5' cgse, (4) 
with 

where j — 14 for AS91 case and 7 = 130 for J86 case. One 
thus has 

Tyi^JgVG ^ ^-^2 ^ ^ eg.Gj-^Rl^^B^,' cos^ an~'-^'. (5) 

According to eq. (Q) , and ((SJ , we have calculated the 
braking indices in different NTVG models, as functions 
of rotation period for typical pulsars with Rq = 1 and 
B12 = 10. The results are shown in Fig.l. It is obvious 
that the braking index could be smaller than 3 for any 
inclination angle in the NTVG model. We find all these 
models can have a braking index as small as 1.4. 

3. Braking indices in emission models where the 
inner and outer gaps coexist 

We simply consider the energy conservation between the 
spindown power and the "dipole radiation -I- inner gap -I- 
outer gap" energy losses, but ignore the interaction be- 
tween these two kinds of gaps. We expect that this simpli- 
fied case will still be useful when more detailed and practi- 
cal consideration is possible in the future. First, however, 
we correct the calculation of braking index^ for the outer 
gap (Cheng Ho & Ruderman 1986; Zhang & Cheng 1997) 
in Paper I, then we present the results calculated for mod- 
els where the inner and outer gaps coexist. 

For a self-sustaining outer gap, which is limited by 
the pair produced by collisions between high-energy 
photons from the gap and soft X-rays resulting from the 
surface heating by the backfiowing primary pairs, the 
potential drop is (Zhang & Cheng 1997) 

A(/.OG - 1.6 X 10i2i?3B-V7r!-io/2i cgse, (6) 

^ Xu & Qiao (2001) made a mistake in calculating the brak- 
ing index for the outer gap, by letting the gap fractional size 
/ = 1 when writing the particle current. 
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and the electric current through the gap 



3.3. SCLF + outer gap 



F = 8.68 X lO^o^Jf r^i^/^i ^ggg^ 



(7) 



thus the energy loss rate due to the outer gap could be 
written as 



(8) 



In the case where the outer gap engine dominates (i.e., the 
inclination angle a = 0, and only the outer gap works), 
Q E'^'^ /Q. ~ fi"^/''; we have therefore a braking index 
n = —5/7 for the outer gap. According to energy conser- 
vation, we obtain the 77 value in a model where inner and 
outer gaps coexist. 



n 



IG+OG 



sin^ a + (5.4 x lQ-^R-^B-^VL-'^/^(t^'^ 
+2.25 X 10^i?,-3 5-12/^17-26/7) cos2 a, 



where Ac^^^ is the potential drop of the inner gap. We 
calculate the braking index, as a function of P, from eq.(^ 
and eq.@. In order to present the correction to Paper I, 
the new braking index calculated for the case where only 
the outer gap works is shown in Fig. 2. 

In the following subsections, we calculate the braking 
indices in the inner-outer gap coexisting models. Three 
kinds of inner gaps are considered, i.e., the normal vacuum 
gap, the near threshold vacuum gap, and the space-charge- 
limited flow. 

3.1. Normal VG gap + outer gap 

Various inner accelerator models have been investigated 
by Zhang, Harding & Muslimov (2000). Here we present 
only the final results of the 77 parameters for the related 
models. Details of the gap potentials can be found in Xu 
& Qiao (2001) and references therein. For the CR-induced 
VG case. 



+2.25 X 10^17-26/7) cos2 a. 

For the ICS-induced VG case, 

T^Jgl+o^ = Sin' a + (1.03 x IQ^^-^^I^ 
+2.25 x 10^0-26/7) cos^ a. 

3.2. NTVG gap + outer gap 

For the CR-induced NTVG case, 



(10) 



(11) 



^NTVG+OG 
'/CR 



sin2a + (1.37/30-"rj-i-76 
+2.25 X 10^17-26/7) cos2 a. 



For the ICS-induced NTVG case, 

^NTVG+OG ^ ^-^2 ^ ^ (6.967-1 X 10217-1' 

+2.25 X 10^17-26/7) cos2 a. 



(12) 



(13) 



For the CR-induced SCLF case, 
^SCLF+OG ^ sin2 a +(38.317-7/4 

+2.25 X 10^17-26/7) cos2 a. 

For the ICS-induced SCLF case, 



SCLF+OG 
"ICS 



sin2 a 



(2.2717-8/13 
2.25 X 10^17-26/7) cos2 a. 



(14) 



(15) 



The corresponding braking indices can be calculated 
from eq.(^ and one of eq.(10)-eq.(15). Fig. 3 shows the 
braking index curves, as functions of pulsar periods. It 
seems that, except for the VG(ICS)+OG model, the en- 
ergy loss of the outer gap dominates. 

4. Conclusion 8c discussion 

The braking index in different NTVG models are calcu- 
lated. The minimum braking indices, ?t-cr^^('^ ~ 0°) = 
1.24 and n[^'s^^(a = 0°) = 1.12, are smaller than 1.4, 
the smallest one observed. This indicates that the NTVG 
model passes the test proposed by Xu & Qiao (2001), 
and that this model may work for radio pulsars. In the 
calculation, typical parameters, e.g., k — pe — Rq — 
B12 = 1 (except B12 — 10 in NTVG models), are em- 
ployed. Nonetheless, we find that the general result does 
not change significantly if these parameters are adjusted 
reasonably. 

It is suggested that the "ICS-VG + Outer" gaps prefer 
to work in the pulsar magnetospheres, since other inner- 
outer gap models imply very small (even a negative) value 
of the braking index, while the observed indices are in the 
range from 1.4 to 2.9. These other models may work if 
much smaller braking indices are observed. For NTVG 
models, there is little difference between the CR-induced 
and the ICS-induced ones, because the energy loss rate 
from NTVG is relatively small compared with that of the 
outer gap. 
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Fig. 1. Set of calculated braking indices, as functions of 

rotation period, for NTVG models. Pulsars arc assumed to 
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here. The curves presented are for pulsars with different 
inclination angles, as indicated in the figure. 
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